Drag forces in classical fields 
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Inclusions, or defects, moving at constant velocity through free classical fields are shown to be 
subject to a drag force which depends on the field dynamics and the coupling of the inclusion to 
the field. The results are used to predict the drag exerted on inclusions, such as proteins, in lipid 
membranes due to their interaction with height and composition fluctuations. The force, measured 
in Monte Carlo simulations, on a point like magnetic field moving through an Ising ferromagnet is 
also well explained by these results. 



In quantum field theory forces between particles are 
induced via their coupling to a quantum field [Ij]. The 
same phenomena arises for fields driven by thermal fluc- 
tuations, for example interactions between inclusions in 
fluctuating membranes |21]. Similarly the Casimir force, 
both quantum and thermal, arises due to the imposition 
of boundary conditions on quantum or thermal fields Q . 
Casimir discovered his force in the quantum context but 
Fisher and de Gennes Q showed that a classical version 
of this effect should be expected for fiuctuating thermal 
fields, such as those for the order parameter of systems 
near a critical point. This critical Casimir effect has only 
recently been measured Q and the technical progress in- 
volved in this experiment may well open up the possibil- 
ity of exploring other aspects of the critical Casimir ef- 
fect, notably dynamical phenomena. In all of the above, 
the effect of the field is manifested by the interaction in- 
duced between two or more particles or surfaces in the 
field. However the presence of the field can also be seen 
by looking at the force exerted on a particle when it is 
not at rest. For instance, electrons moving in materi- 
als induce a local polarization known as the polaron fil 
which modifies their dynamics. A frictional Casimir force 
is also induced by the uniform motion of a conductor in 
a volume of blackbody radiation which is in equilibrium 
in the rest frame of a cavity containing the radiation 7] . 

In this letter we show that for classical fields, in a lab- 
oratory rest frame, a drag force is present on inclusions 
linearly coupled to the field, and which move at constant 
velocity. The underlying physics is caused by a polaron- 
like phenomena (see Fig. ([T])) which we generalize to a 
range of statistical field models arising in the study of soft 
condensed matter systems. A key point in this analysis 
is that we examine the effect of the dynamical models 
commonly used to study soft condensed matter systems 
on the drag forces on inclusions. 

As a concrete example of the class of problems we will 
address we start by studying drag forces in the Ising fer- 
romagnet model on a 0? dimensional cubic lattice with 
Hamiltonian 
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FIG. 1: Contour plot of the magnetization profile (polaron) 
for the 2d Ising model about a local magnetic field at a single 
point moving with velocity v. (high temperature phase: /3 = 
1, J = 0.4, h = 6.66). 



neighbor spins on a square lattice of spacing a, and h 
a magnetic field at the position of an inclusion iq which 
moves in the direction z with a velocity f , so io{z) = vt/a 
(where we take the integer part). The system evolves 
in the following manner, the underlying unit of time is 
a Monte Carlo sweep where N (the system size) ran- 
domly chosen spins are examined and flipped or not ac- 
cording to the dynamical rules used: (i) dynamics not 
conserving the total magnetization - Glauber dynamics 
- a single spin is chosen and is flipped with probability 
Pf = -I- exp(^AiJ)], where /3 is the inverse tempera- 
ture and AH the energy change associated with the spin 
flip (ii) a form of Kawasaki dynamics conserving the to- 
tal magnetization, here two spins of opposing sign are 
randomly selected and flipped with the probability pj. 
After T ^ a/v units of time the inclusion is moved one 
step in the z direction. The instantaneous force is defined 



Here J > is a ferromagnetic coupling between nearest as f{t) = —AE{t)/a — —h{Sig+i{t) — Si„-i{t)) /2a (the 
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energy change Ai? being computed from the symmetric 
discrete derivative on the lattice). 

Shown in Figs. ([2^) and ^jp) are the steady state val- 
ues obtained for the average value of the force f{v), for 
Glauber and Kawasaki dynamics respectively, as a func- 
tion of the velocity v for the one dimensional Ising model 
with parameters, P — 1, J — 1 and ft, = 10. Similar re- 
sults are found in two dimensions. All numerical results 
have similar characteristics, the force is always in the op- 
posite direction to the direction of movement, so it is a 
friction, the force is linear for small v and for large v the 
force decays as 1/v. 

Fig. dl]) is a plot of the average magnetization around 
the inclusion seen from its rest frame in the steady state 
regime with Glauber dynamics. For zero velocity we see 
that the magnetization profile is spherically symmetric 
and there is no net force, but as the velocity is increased 
the magnetization profile is distorted and takes a differ- 
ent form ahead of and behind the inclusion. As v in- 
creases the smearing out of the magnetization profile is 
increased, but the amplitude of the deformation dimin- 
ishes. These results indicate that it is the magnitude of 
the deformation of the local field by the inclusion along 
with the asymmetry induced by the particle motion that 
induces the force. Once the velocity becomes too large 
the system does not have time to react to the presence of 
the inclusion and the frictional force weakens, this shear 
thinning like cross over should be quite generic. These 
results are qualitatively similar to those of [8^ for the in- 
teraction of an MFM tip passing over a magnetic sample, 
with dipolar interactions between the tip and the sam- 
ple, and with spin dynamics described by the Landau- 
Lifschitz-Gilbert equation. 

To understand these results we consider the general 
over-damped dynamics for a scalar field (j). We denote 
positions in space by r = (x, z) where the coordinate z 
is in the direction of the motion. We write the Hamil- 
tonian density H as two parts, the first a bulk term 
'Ho[0] and the second the energy due to the inclusion 
S{x.)5{z — vt)Htr['P]- We thus assume that the interaction 
between the inclusion and the field is localized about the 
point (0, vt). The energy of the system is thus 
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H = / {Ho[cl)] + S{x)S{z - Z)ntr[<P]) dr, 
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where Z = vt is the inclusion's position in the direction z. 
The instantaneous force acting on the inclusion is simply 
given by 



/ = - 



dH 
'dZ 



dHtrimz)) 

dZ 
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and when the derivative is not continuous we implicitly 
take its average value about the particle position (cor- 
responding to the symmetric discrete derivative used in 
the Ising model simulations). We now take a general 
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FIG. 2: Dashed lines: average drag force / in the Id Ising 
model as a function of v for Glauber (a) and Kawasaki dy- 
namics (b). Solid lines are the fits of model A (a) and model 
B (b) dynamics for the Gaussian ferromagnet. 



over-damped dynamics for the evolution of the field 
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where i? is a dynamical operator defining the under- 
lying over-damped dynamics. For instance i?(r, r') = 
i?(r— r') = (5(r— r') for model A dynamics (non-conserved 
order parameter) and R{r — r') = — V^(5(r — r') for model 
B dynamics (conserved order parameter). The noise 
term is taken to be Gaussian white noise in time and 
with spatial correlations that respect detailed balance, 
i.e. (r;(r, t)r/(r', t')) = 2TS{t - t')R{v - v'). 

Consider the case of free scalar fields with a linear cou- 
pling between the inclusion and the field i.e. 
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Ho = -0(r)A0(r), Ut 



-K(j){r) 



(5) 



where both A and K are linear operators which, for the 
purposes and applications here, we will assume to be 
self adjoint. Working in the coordinate system comov- 
ing with the inclusion, the evolution equation for average 
value of the field takes the form 



9(/)(r) (90(r) 
dt dz 



-RA(l){r) + RK{r). 



(6) 
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In the steady state we can neglect the time derivative 
in the above and solve the resulting equation by Fourier 
transforming Q (where g(k) — J dxexp(— ik-x)g(x)) to 
obtain 



<^(k) = ^ 



R{k)K{k) 



A(k)i?(k) - ik^ 
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A key point of Eq. ^ is that when v the average 
magnetization depends on the dynamical operator R, it 
is only when v = that (p takes its equilibrium value. 
Using this we find the average force to be 
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A(k)i?(k) - ik^ 
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where d denotes the spatial dimension. The force thus 
depends on the particular form of the dynamics via R as 
well as the static quantities A and K. For small velocities 
we can write this result in terms of an effective friction 
coefficient f = —Xv where 
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{2TTYdJ A^k)R{k) 
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with k = \k\ and we have assumed that the operators 
A, K, and R are rotationally as well as translationally 
invariant. For large velocities we find that 



dkR{k)k'^{k). 
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The limiting forms for the average force for small and 
large velocities thus take a universal form as long as 
the integrals in Eqs. (jU and ((T0| exist. These in- 
tegrals can be regularized at large k by an ultraviolet 
cut off. However if the insertion has a finite size, then 
we can approximate this by replacing the delta func- 
tion in the above computation by a Gaussian profile of 
width a corresponding to the insertion size. This gives 

K{k^ a) — K{k, 0) exp ^— -^-j-^ , but this simply sets an- 
other cut-off fcmax ~ 27r/a. 

To see what the generic form of the induced drag is 
we introduce a time scale t for the dynamics writing, 
on dimensional grounds, i?A = a^^^ k^^^ /t, and where 
for simplicity (however see the discussion later) we as- 
sume that flQ = a is the cutoff of the field theory as well 
as the inclusion size. We also write A = /j,ag~^fc* and 
K = ^'a^k". A critical dimension dc — 2S + p — 2a — 2 
is found such that if d > dc, and we assume that the 
friction is dominated by the ultra-violet divergence, we 
find A ^ ^t/gq. However if d < dc, but there is a mass 
term m wfiich regularizes the theory at the infrared level 
then we find A ^t{£,/ a^Y""''- / a^, where ^ = 1/m is 
the correlation length of the field. Hence when d < dc 
the typical friction can be expected to be orders of mag- 
nitude greater than that for systems such that d > dc 



if the correlation length is significantly larger than the 
inclusion size. When d < dc, but the system is critical 
i.e. m — (corresponding to the Gaussian approxima- 
tion for continuous phase transitions such as the para- 
ferromagnetic transition or critical demixing in lipid bi- 
layers) the integral is infrared regularized at kmin 2tt/L 
where L is the linear dimension of the system and we find 
A ^ l^r{L/aor^-'^/al 

Note that this critical behavior can even be present 
when the free field theory is not critical (m ^0); in this 
case dc = p — 2q; — 2. This strange situation can occur 
when the long distance modes of the theory are weakly 
damped (p large for A; <C 1) or when the field coupling 
to the inclusion position is long range (non local a < 0). 
Critical behavior can also be seen in the expression Eq. 
(|S]) at small v. li S + p > 1, which is the case for the 
examples given in this letter, and d < dc we find that 

7(u) ^ x;l+P+a-i . 

Gaussian Approximation to the Ising Ferromagnet: to 
understand the numerical results on the Ising model, we 
take, in the disordered phase, the Gaussian model A — 
k^+m? . The coupling to the inclusion is linear so K — h. 
For Glauber dynamics we take the model A dynamics 
R — 1 and for Kawasaki dynamics we take model B with 
R — k^. When m = 0, dc — 2 for model A and dc = ^ for 
model B. Shown in Figs ((T^) and ((TJj) are the simulation 
results for the Id Ising model with Glauber and Kawasaki 
dynamics fitted with the formula Eq. ([8]) using model A 
and B dynamics for the Gaussian model, we see that the 
fits are excellent in both cases. 

Drag in lipid membranes: SafFmann and Delbriick (SD) 
[loj computed the diffusion coefficient for cylindrical in- 
clusions in membranes where the membrane is treated 
as a continuous medium. This was done by computing 
the drag by taking into account the coupling of the 2d 
hydrodynamic field of the membrane to that of the exter- 
nal medium using the nonlinear Navier Stokes equations. 
The diffusion constant is then obtained via the Stokes 
Einstein relation and is found to depend on the protein 
radius a as I? ^ ln(a). However a recent experimental in- 
vestigation [Ol found that the diffusion constant behaves 
as D ^ 1/a. Several factors can change the dependence 
of the diffusion constant on a, for instance with the same 
hydrodynamical model frictional coupling of the mem- 
brane to a solid surface leads to a behavior D ~ 1/a^ 
[l^ . Deviations from the SD formula have also been at- 
tributed to modifications of membrane thickness, due to 
hydrophobic mismatc h. Il3l| and also the modification of 
membrane curvature |14| caused by the protein. 

In the spirit of these latter approaches, first consider 
the case where the parameter corresponds to the height 
of the membrane. In the Monge gauge, for small fluctu- 
ations, A(A;) ~ Kfc^ + crfc^ where n i s the bending rigidity 



and a is the membrane tension 15|. A simple model for 



the coupling between a protein and a membrane is a lin- 
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ear coupling between the membrane curvature and the 
protein (a protein favors a positive of negative curvature 
of the membrane), this means that K = hk?. The dy- 
namics of membrane height fluctuations are dominated 
by hydrodynamic interactions which can be shown to lead 
to a nonlocal dynamical operator R = l/Arjk, where rj is 
the viscosity of the surrounding medium 16]. 

To model drag in membranes we need to know how the 
magnitude of the coupling K depends on the size of the 
insertion a. To find the dependence of if on a we assume 
that the energy induced by the inclusion has the form of 
a line tension for small a, that is we impose 



2(27r)'' 



A{k) 



27ra7, 



(11) 



where 7 is the line tension. We set K = k'^h and assume 
that the integral is dominated either by an ultra-violet 
divergence, or when this is not present by an infrared 
divergence regulated by fc,„i„ = rn = 1/^, yielding, 
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d+2a~S+l 



for d - 
for d 



2a~S>0 
-2a-S <0 



a/ In 



for d+2a 



0. 



(12) 
(13) 
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Using this result to compute the friction for curvature 
coupling to membrane height fluctuations we find X ^ a^. 
This means that the frictional force is relatively weak, but 
if it were dominant one would estimate, via the Stokes 
Einstein relation, that the diffusion constant would scale 
as D = ksT/X ~ 1/a^. Protein coupling to local lipid 
composition (either chemical or order) can induce inter- 
actions between proteins 17[. Let us take (p to be ei- 



ther the local order parameter specifying the local phase 
(solid, liquid, gel) for single lipid membranes or the lo- 
cal lipid composition for membranes composed of lipid 
mixtures 18[, and assume it couples linearly to the pro- 
tein position, i.e. K = h \s constant. The experiments 



of [ll| on the diffusion constant of small peptides were 
carried out in lipid membranes composed essentially of 
a single lipid type SOPC, the order parameter coupled 
to the protein has no reason to be conserved and so we 
assume dynamics of model A type. Applying our compu- 
tations for the magnitude of h along with the dynamical 
result for the friction Eq. © we find A ~ a. Thus for 
small insertions this drag force is expected to be larger 
that the height fluctuation induced force and if this drag 
dominates we obtain the behavior D ^ 1/a - the scaling 



found in 11|. It is thus possible that, for small membrane 
inclusions, the dominant drag force is due to interaction 
with an order parameter of the surrounding lipids. 

We have shown how inclusions moving at constant ve- 
locity in classical fields experience a drag force. This drag 
force is generated by a polaron-like perturbation of the 
surrounding field, at large velocity the polaron is weak- 
ened as the field cannot react sufficiently quickly to the 



inclusion. The results given here are for free fields, but a 
rich phenomenology arises from a relatively straightfor- 
ward analysis. The wealth of the theory comes in part 
from the range of field theories that describe soft-matter 
systems and from the range of dynamical models that 
are present. There are clearly a wide range of directions 
for further study, including (i) the effect of interacting 
fields, (ii) extensions to more realistic dynamics and (iii) 
the possibility of observing the drag behavior predicted 
here in experiments. For example this phenomena could 
possibly be studied using optical tweezers to pull an in- 
clusion, such as a colloid, through a binary fluid mixture, 
in particular at its critical point. It would also be inter- 
esting to see how protein coupling to local membrane 
composition affects the diffusion constant for proteins in 
multicomponent membranes. 
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